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SUPER DIRECTIONAL BEAMFORMING DESIGN AND IMPLEMENTATION 

Cross-reference to Related Applications 

This application is related to application U.S. Serial 

No.: , by Andrea et al., filed on October 22, 1999 and entitled 

5 System and Method for Adaptive Interference Canceling (Attorney Docket No. 670025- 
7010), hereby incorporated by reference. 
Field of the Invention 

This invention relates to signal processing, and more particularly, to 
processing the signals received by an array of sensors in order to minimize the amount 
1 0 of noise received by the array when the array is being used to receive a desired signal. 

Background Of The Invention 

Beamforming is a term used to designate the operations associated with 
forming spatial sensitivity pattern for an array of sensors. Classical beamforming is 
defined as "delay and sum beamforming". In delay and sum beamforming, a source 

15 transmits a wave that propagates and arrives at an array of sensors at different times, 

depending on the source direction and the array geometry. The outputs of the sensors of 
the array are delayed, to compensate for the delay in time of arrival of the source's 
wave, which originated from the preferred direction, and summed, to provide a classical 
directional beamformer output. The effect of sources that are located at directions other 

20 than the preferred direction (referred to as the looking direction) is reduced by the 
beamforming process, resulting in maximum sensitivity of the process towards the 
preferred direction. 
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The array of sensors can be, for example, an array of microphones 
receiving an acoustic sound source. The beamforming process can be used to map 
sound sources (in a sonar system for example), or to emphasize a sound source whose 
direction is known, by modifying the compensating delays and "steering" the look 

5 direction of the array. The beam-width — usually defined as the difference between the 
two angles, in which the output energy is reduced by 3dB relative to the beam center — 
depends on the array length, frequency of the received signal and propagation speed of 
the received signal (in our example the speed of sound). For many practical purposes 
the beam-width of the array will not be sufficiently narrow, and enlarging the array 

1 0 length is not desired. For those cases a more directional beamforming process is 
required. 

Moreover, while delay and sum beamforming, does not provide optimum 
noise reduction. If the sensors' outputs are filtered (a different filter to each sensor) and 
the outputs of the filters summed, one can obtain a different shape of the beamfonner 

1 5 output and improve noise reduction. With a careful design it is possible, for example, to 
create a null (zero reception) towards a given direction. If a noise source's direction is 
known and a null is placed in that direction, improved noise reduction can be realized as 
compared to the noise reduction of the classic delay and sum beamfonner. 

Two basic approaches have been developed to obtain optimum 

20 performance of a beamfonner in the presence of noise. The first one, presented in 

Monzingo and Miller - Introduction to Adaptive Arrays (Wiley, New York) pp. 89-105 
and 155-216 shows that if a filter is created for each sensor that for each frequency will 
have gain weights of 
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C 'v 
= r~ (1) 



10 



the output of the beamformer will have optimum performance in terms of noise 
reduction. The above weights will maintain a unity gain at the look direction (no 
distortion of the desired signal) while providing minimum energy at the output. The two 
assumptions (minimum energy and no signal degradation) will result minimum noise at 
the output. In Eq. (1) C is the noise covariance matrix and it may be expressed as: 



C = £{y*y T ) 



(2) 



where 



y T =[y»(0 y 2 (f)..yn(i)l (3) 

is the noise measurement at the elements, and v is the steering vector towards the look 
direction, expressed as: 



15 



v = 



(4) 



where 



T0-x(n-l) are the steering delays introduced to elements 0-n respectively by a target 
originated at the look direction. Further, the filtered elements approach was extended by 
20 Frost (O. L. Frost, III, "An Algorithm for Linearly Constrained Adaptive Array 
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Processing" Proc. IEEE, vol. 60, no. 8, pp. 926-935, Aug. 1972.) to provide an 
adaptive beamformer in which the weights would adapt themselves so that they 
converge to provide the optimum solution. 

The second basic approach to obtain optimum beamformer performance 
5 was developed by Griffiths (L. J. Griffiths and C. W. Jim, "An Alternative Approach 
to Linearly Constrained Adaptive Beamforming," IEEE Trans. Antennas PropagaL, 
vol. AP-30, no. 1, pp. 27-34, Jan. 1982.) who suggested using a Noise Canceling (NC) 
approach to the optimum beamformer problem. In his approach the adaptive coefficient 
are updated by the Least Mean Squares (LMS) algorithm. Griffiths proposed using the 

10 elements' signals to obtain a main channel, in which both the signal and the noise are 
present, and reference channels, in which only noise is present (i.e. which are signal 
free). The main channel can be generated through one of the elements alone, or through 
classic delay and sum beamforming. The reference channels can be generated through 
the subtraction of one element from another, or by forming any other linear combination 

15 of elements that would provide a zero output at the look direction (i.e. the signal 
direction). The main channel and the reference channels are utilized by an adaptive 
LMS Widrow filter to obtain an optimum beamformer (see Adaptive Noise Canceling: 
Principals and Applications - Widrow, Glover, McCool - Proc. IEEE vo!63 no.12 
1692-1716, Dec 1975). In this adaptive beamformer each reference channel is filtered 

20 (i.e. each channel signal is convolved with a set of filter coefficients), the filtered 

channels are summed together to obtain the noise estimation, and the noise estimation is 
subtracted from the main channel to provide a noise free signal. The filter coefficients 
in the Griffiths solution converge to 
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w op t=Cr , p (5) where 

C is the noise covariance matrix and p is the correlation vector between the beam 
output and the reference channels. Note that with this approach the steering is done 
through the creation of the reference channels and the beam, so there is no steering 
5 vector towards the look direction in equation (5). Griffiths showed that, for an n 
elements system, if one creates n-1 reference channels, the LMS approach would 
-converge to the same optimum solution as Frost.- — 
Objects and Summary of the Invention 

It has been recognized that while the two approaches to optimum 
1 0 beamforming discussed above were primarily developed to provide an adaptive 
solution, they also teach us what the optimum solution would be given the noise 
covariance matrix. A non-adaptive approach, in which predetermined filters are 
designed and used, is sometimes more appealing than an adaptive approach. The fixed 
beam (non-adaptive) approach is much less computationally intensive, it is much less 
15 sensitive to leakage of the desired signal to the reference channels and it does not give 
rise to distortion in the desired signal. Also, the fixed approach has the potential to 
handle some types of noises better than an adaptive process, such as reverberation and 
diffused low noises. On the other hand, one may not want to give up the adaptive 
process, because it provides the best immunity to significant directional noises. A 
20 hybrid system that uses both adaptive and non-adaptive techniques provides a system 
which realizes the advantages of both techniques. 

Further, it has been recognized that while the above described optimum 
beamforming techniques provide the solution given the noise covariance matrix, they do 
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not show how to determine this matrix for a particular noise scenario. Also, the 
equations show how the required weights for each frequency can be computed, but they 
do not show how to implement the time domain filters that will approximate the 
weighting solution. The prior work in this area does not discuss how such time domain 
5 filters would be designed or implemented in a combined adaptive/non-adaptive 

beamforming system. Moreover, there is no teaching as to techniques for overcoming 
differences in the elements' sensitivity, phase, or the influence of packaging and other 
mechanical interferences on the performance of the fixed beam. 

In view of the above considerations, it is an object of the invention to 
10 provide a sensor array beamforming system capable of optimal noise reduction 
performance. 

It is another object of the invention to provide a simple and easy method 
to design optimal filters in a sensor array beamforming system. 

It is still another object the invention to provide a simple and easy way to 
1 5 implement the optimal system in as a fixed solution system or as a combined fixed and 
adaptive system. 

It is yet another object of the invention to provide a method to design 
optimum filters for a sensor array beamforming system that would take into 
consideration the specific characteristics of the sensors (microphones for example), and 
20 other mechanical or acoustical features that influence the performance of the array. 

In order to realize the above objects of the invention and overcome the 
drawbacks of prior systems, the invention provides a sensor array receiving system 
which incorporates one or more filters that are capable of adaptive and/or fixed 
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operation. The filters are defined by a multiple of coefficients and the coefficients are 
set so as to maximize the signal to noise ratio of the receiving array's output. In one 
preferred embodiment, the filter coefficients are adaptively determined and are faded 
into a predetermined group of fixed values upon the occurrence of a specified event. 
5 Thereby, allowing the sensor array to operate in both the adaptive and fixed modes, and 
providing the array with the ability to employ the mode most favorable for a given 
operating environment. In another preferred embodiment, the filter coefficients are set 
to a fixed group of values which are determined to be optimal for a predefined noise 
environment. 

10 Thus, reference is made to application U.S. Serial 

No.: , by Andrea et aL, filed on October 22, 1999 and entitled System 

and Method for Adaptive Interference Canceling (Attorney Docket No. 670025-7010), 
which, together with the documents and patents and patent applications cited therein are 
hereby incorporated by reference; the present invention may be used in conjunction 

15 with embodiments disclosed or discussed in Andrea et al. and/or in the documents, 

patents and patent applications cited in Andrea et al. (and incorporated herein), e.g., the 
"Superbeam" technology of this invention can be used in conjunction with "DSD A" 
technology in embodiments disclosed or discussed in Andrea et al. and/or in the 
documents, patents and patent applications cited in Andrea et al. (and incorporated 

20 herein). 

Brief Description of the Drawings 
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The following detailed description, given by way of example and not 
intended to limit the present invention solely thereto, will best be appreciated in 
conjunction with the accompanying drawings, wherein like reference numerals denote 
like elements and parts, in which: 

5 Fig. 1 is block diagram of a filtered input type beamforming system in 

accordance with the present invention. 

Fig. 2 is a block diagram of a filtered references type beamforming 
system in accordance the present invention. 

Fig. 2A is a flowchart which shows an illustrative procedure for 
10 designing and implementing the fixed filtered references approach. 

Fig. 3 is a flowchart showing an illustrative procedure for generating 
fixed filter coefficients through the use of simulated noise and an actual adaptive system 
positioned in an an-echoic chamber. 

Fig. 4 is a flowchart showing an illustrative procedure for generating 
15 fixed filter coefficients through the use of simulated noise, a microphone array 

positioned in an an-echoic chamber and an actual adaptive system positioned outside an 
an-echoic chamber. 

Fig. 5 is a flowchart showing an illustrative procedure for generating 
simulated noise and using the simulated noise to generate fixed beamformer 
20 coefficients. 

Detailed Description 

The following description will be divided into four parts. Part one will 
detail a method for designing and implementing fixed beam optimal filters based on the 
filtered input approach. Part two will detail a method for designing and implementing 
25 fixed beam optimal filters based on the filtered references approach. Part three will 

detail a hybrid system that includes both a fixed solution and an adaptive one. Part four 
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will detail two alternative approaches to the design and implementation of fixed beam 
filters. 



1 - A metho d of designing and implementing fixed beam optimal filters based on the 
5 filtered input approach. 

Fig. 1 is block diagram of a filtered input type beamforming system in 
accordance with the present invention. As can be seen from the Fig. 1, N microphones 
IOi.n are conditioned and sampled by signal conditioners 12^. The microphones' 
samples are respectively stored in time tapped delay lines 14|. N and filtered by filters 
10 16i- N via convolvers 1 8i. N . The output of the filters is summed up via an adder 20 to 
provide a fixed beamformer solution. In the Fig. 1 embodiment, the invention provides 
a method for creating the noise covariance matrix and then using equation (1) to 
actually design the coefficients for filters 16t= N . It should be noted that the equation 
provides us with the required coefficients in the frequency domain. The time domain 
15 coefficients are obtained from the frequency domain coefficients. 

To determine the optimal solution for any given scenario we first define 
the scenario in terms of the spatial distribution of the interfering sources (directions and 
relative intensity). For each of the interfering sources we assume a far field model. Let 
us consider an array of M identical omni-directional sensors with a known arbitrary 
20 geometry measuring the wave-field generated by a single far-field source. Let denote 
the location of the f-th sensor, where n =[Xj, y t , z{] and let ^ and 0 denote the azimuth 
and elevation angles of the radiating source, respectively. 
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Let us now define a differential delay vector, which expresses the delay 
in time of arrival of the interference wave front to the various elements: 

where the first sensor serves as a reference and the delays are measured relative to it. 
5 The signal "Direction Of Arrival" vector for the far field case is given by: 





V 




sin(#)cos(j0) 


k = 


*, 




sin(#)sin(^) 








cos(0) 



(6) 



10 



Let us define a distance matrix R between the sensors of the array 



»2-r, 



(7) 



The time delay between any two sensors is equal to the projection of the distance vector 
between them along the k vector divided by the wave propagation velocity (sound 
1 5 velocity for example). Consequently, the delay vector can be expressed as follows: 

Rk 



T - — 



(8) 
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where c is the wave velocity and the matrix R is composed of the distance vectors 
between all the sensors and the reference sensor. More explicitly, for sensor j we can 
write 

5 r„. =[x iy cos(<9)sin( ( z>) + >' iy sin((9)sin(^) + z i; cos(^]/c (9) 

Assuming that interference i has an amplitude of Si and a Direction Of Arrival vector of 
kj then its measurement by the array can be expressed as the source steering vector 
multiplied by the source amplitude 

10 

The contribution of source i to the noise covariance matrix is expresses as: 

15 C,.=£{y,y/} (11) 

Since bj is deterministic and we assume stationary sources where s 2 is the power of the 
source i the above equation is reduced to 

c,=y f y/ Ola) 

20 Under the assumption that the interferences are uncorrelated we can write 
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12 

C = £C t (12) 

If we assume that there is an additive uncorrelated noise (spatially distributed white 
noise) n to each of the sensors we obtain 

C = «I + £C. (13) where 



I is the unity matrix with a size of [MxM]. 

So far we obtained the noise covariance matrix for a predetermined noise 

1 0 environment. In order to use equation (1 ) we need to calculate the steering vector v. 
This steering vector expresses the look of the array towards a defined direction. The 
steering vector v is the conjugate of the vector already expressed in equation (4) and it 
is calculated in the same way as the steering vector of the noise sources (see (8) and (9)) 
where <f> and 6 are the azimuth and elevation of the look direction, respectively. 

1 5 It should be noted that a far field model for the noises was used to obtain 

the above equations. It is not necessarily desirable to use a far field model for the target 
(desired signal). For example, one may want to implement a focusing effect on the 
target in near field situations. Such an effect can be obtained by manipulating the 
steering vector accordingly. 

20 The fixed solution technique of Fig. 1, using equation (1), provides a 

way to calculate the gain weights of each sensor in an array for each frequency. More 
specifically, for each frequency of interest the system of Fig. 1, equation (1) is solved to 
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yield one weight for each filter (w op t is a vector with the number of elements being 
equal to the number of sensors). Thus, if it is desired to obtain the optimum weights for 
ten frequencies, for example, equation (1) is solved for ten frequencies and each filter 
16i-n is then defined by ten frequency domain weights - the set of frequency domain 

5 weights for each filter defining the filter's frequency domain response. 

Once the frequency response for each filter is determined, it is necessary 
to design the time domain filters to provide the determined frequency response. If the 
weights (or "gains") are real numbers - meaning that the desired filter has a linear phase 
. we can use the weights with any of the well-known methods to design the filter for 

10 each sensor. For example, a Remez Exchange Method can be used. For simple cases 
such as when the array is linear and the noise sources are positioned in a symmetric 
structure around the look direction, the gain weights would be real numbers. If the gain 
weights are complex numbers, such as when the noise structure is not symmetric, the 
required filter will not have a linear phase. For these cases one can feed the weights for 

15 each filter to an IFFT (Inverse Fast Fourier Transform) procedure to obtain the time 
domain function that would provide the desired frequency response and phases for the 
filter. 

2 - A method for designing and implementing fixed beam optim al filters based on the 
20 filtered references approach . 

Fig. 2 is a block diagram of a filtered references type beamforming 
system in accordance the present invention. As can be seen from Fig. 2, N microphones 
26i_ N are conditioned and sampled by signal conditioners 28i- N . The microphone outputs 
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are processed by a delay and sum beamformer 30 to provide a beam channel, and by a 
reference channel processor 32 which is typical of an LMS beamforming system. As 
shown, the beam channel may be formed via the classic delay and sum beamforming 
process on the inputs, however the alternatives include any linear combination of sensor 
5 outputs that will provide a maximum towards the looking (listening) direction. The 
reference channels are processed such that a null is placed towards the looking 
direction. It may be obtained by subtracting one microphone form the other, or by 
forming some other linear combination of sensor outputs. The output of the reference 
channels is respectively stored in tapped delay lines 34i. L (L may or may not be equal to 

10 N) and filtered by filters 36|. L via convolvers 38i- N The filtered reference channel 
output is summed via an adder 40 and subtracted via a subtractor 42 from the 
beamformer output as delayed by a delay line 44. This structure is typical to adaptive 
beamformers, where the reference channels are filtered by adaptive filters and then 
summed and subtracted from the delayed main beam signal. In our case, the filters are 

15 fixed (non adaptive) and pre-designed. The method is highly practical in systems that 
already have the structure of an adaptive beamformer, which can be applied to both the 
adaptive solution and the fixed solution. 

In the filtered references embodiment of Fig. 2, the filters' coefficients 
are designed and determined using equation (5). More particularly, the noise covariance 

20 matrix is determined and then used in equation (5) to determine the filter coefficients. 
As was the case in the filtered input embodiment of Fig. 1, equation (5) provides filter 
coefficients in the frequency domain and it is necessary to obtain the time domain 
coefficients from the frequency domain coefficients. 
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Equation (5) is expressed as 

w 0/ ,,=C _I p where 

C is the noise covariance matrix as measured by the reference channels, and p is the 
correlation vector between the main channel (beam) output and the reference channels. 
5 We obtain the noise covariance matrix using techniques that are similar to those used in 
the filtered inputs approach. The difference is that we need to obtain the noise received 
as it appears in the reference channels, and hot as it appears at each sensor. To do this, 
we first obtain the contribution of each noise source to each sensor (the same y\ that we 
obtained in the previous method), and then find the contribution of each noise source to 

10 each reference channels. The reference channels are generally relatively flat sensitivity 
patterns having nulls pointing to the array look direction. The reference channels are 
created using linear combinations of the elements* outputs after they have been steered 
to the look direction. For example xi+x 2 -(x 3 +X4) may be a reference channel after the 
inputs (denoted as x n ) have been appropriately delayed to compensate for the look 

15 direction. These relations can be expressed as a nulling matrix N (note again that in 
order to guarantee a signal free reference the sum of the elements of each row in the 
matrix should be 0). 

Example for nulling matrix for an array of four microphones and three reference 
20 channels is 



H-I 
4 



1 1 -1 -1 
■3 11 1 
111-3 



(14) 
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Note also that for an n elements array only n-1 independent nulls can be created. If we 
denote v as the steering vector to the look direction than we can obtain the contribution 
i of the a noise source i to the reference channels through the following equation: 

x, =fidiag(Y)y t (15) 



where diag(v) is the diagonal matrix which elements are the element of the vector v (for 
broad side array diag(v)=I - the unity matrix), y, is the interference contribution of 
10 noise source i measured by the array elements as described above, N is the a Nulling 
matrix used to create the reference channels and x f is the contribution of interference of 
noise source i as measured by the reference channels. Through equation (15) the 
contribution of a noise source is "transferred" from the array elements to the reference 
channels. 



15 



The overall noise measured the reference channels is the sum of the 
noise contributed by each interference. 



x = £x, (16) 

20 

where x is the noise measured at the reference channels. 

The contribution of each x, to the noise covariance matrix is expressed as 
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c ( . 0 7 ) 

As in the case of equation (11), since x is a multiplication of a stationery signal by a 
deterministic one (the steering elements) the equation is reduced to 
5 C,. = x, X ; (17a) 

Under the assumption that the interferences are uncorrelated we can write 

10 If we assume that there is an additive uncorrelated noise n (spatially distributed white 
noise) to each of the sensors we obtain 

c=«i+£c f < 19 ) where 

15 I is the unity matrix with a size of [M x M]. 

We now need to find the correlation vector p. This vector expresses the 
correlation between the beam signal and the reference channels. The correlation vector 
p is given by: 

20 P = I> ( 20 ) 

where pi is given by 

p,. =beam.x. (21) 
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and beam i = v'y, (22) 

After obtaining both C and p equation (5) is used to find the gain weights for each 
frequency. The practicality of obtaining the weight for a series of discrete frequencies 

5 and the actual design of the filters was demonstrated above in relation to the filtered 
inputs method of Fig. 1. 

An illustrative procedure for designing and implementing the fixed 
filtered references approach is shown in Fig. 2A. As can be seen from the figure, the 
first steps are to define the desired noise scenario, the array configuration and frequency 

1 0 range and resolution (step 50), and to initialize certain variables to be used in the 

procedure(step52). Next, the contribution of a first noise source to the noise covariance 
matrix - at the array output - is computed (step 54). The noise source's contribution to 
reference channel noise covariance matrix is then computed on the basis of the source's 
contribution at the array output, the nulling matrix and the steering vector toward the 

1 5 array look direction (step 56) is computed, and the correlation vector between the beam 
signal and the reference channels for the source is determined (step 58). 

At this point a determination is made as to whether each source has been 
considered in steps 52-58 (step 60). If not all noise sources have been considered, a 
count variable is incremented (step 62) and steps 52-58 are performed for the next noise 

20 source. If all noise sources have been considered, the contributions of each noise source 
to each reference channel are summed to generate a reference channel covariance 
matrix and the beam/reference channels correlation vectors are added to determine a 
beam/reference channel correlation matrix (step 64). Once the reference channel noise 
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covariance matrix and correlation matrix are determined for a particular frequency 
under consideration, a filter coefficient corresponding to that frequency is determined 
for each channel according to equation (5) (step 66). 

A determination is then made if each desired frequency has been 
5 considered (step 68) If not all frequencies have been considered, a count variable is 
incremented (step 70) and steps 52-68 are performed for the next frequency. If all 
frequencies have been considered, a filter design program is used to obtain the filter 
time domain coefficients that approximate the desired response as defined by the 
frequency domain coefficients determined in step 66 (step 72). 

10 

V A hybrid system that includes both a fixed solution a nd an adaptive one. 

Adaptive systems are designed to provide the optimum solution to the 
noise environment at any time. Using the reference channel type approach, an adaptive 
system measures and studies the noise sources through the reference channels and 

1 5 subtracts it utilizing LMS filters. A major problem of an adaptive system is the leakage 
problem. The desired signal "leaks" into the reference channel nulls due to differences 
in the sensors' sensitivity and phases, or due to mechanical imperfections. The leakage 
of the desired signal into the nulls causes the system to try and cancel the desired signal 
as though it was noise, and thereby causes distortion in reception of the desired signal. 

20 One way to prevent signal distortion due to leakage is by blocking (or freezing) the 

adaptive process when a strong desired signal is detected, and thus prevent the adaptive 
process from attempting to cancel the desired signal. However, regardless of the logic 
of the adaptive process blockage, blocking has the effect of locking the noise reduction 
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filters on the solution existing immediately before blockage commenced, resulting in 
the filters losing their relevancy in time. 

In order to overcome the ameliorate the problems associated with 
leakage and blocking, the present invention provides a system in which the filters' 

5 coefficients drift form their adaptive solution into a pre-designed fixed solution. The 
system initializes its filters' coefficients with the fixed pre-designed solution and fades 
into the fixed solution whenever the adaptive process is blocked. The drifting 
mechanism is implemented in the following way: let Wj(n) be the i-th coefficient of an 
adaptive filter at time n, and let w(0) be the fixed value of that filter coefficient, then 

10 Wi(n+l)=wi(n)*y+w,(0)*(l-y) (23) where 

y determines how fast the filter will converge into its fixed solution. 

The drifting process of the invention serves another purpose. It has been 
1 5 shown that the adaptive process may explode (or diverge) due to numerical problems 
when the process is performed by a fixed-point processor (see Limited-Precision 
Effects in Adaptive Filtering - John M. Cioffi - IEEE Transactions on Circuits and 
Systems vol cas-34 no. 7, Jul 1987). To prevent such a divergent breakdown, it is 
sometimes useful to apply a "leaky filter". A leaky filter multiplies its coefficients by a 
20 number smaller than one before they get updated, thus preventing divergence due to 

numerical problems. Although the leaky process does not allow the filter to converge to 
the optimum solution, it prevents mathematical divergence. 
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The use of the decaying process proposed here will eliminate the need to 
use a strong leaky process (or any leaky process) since whenever the adaptive process is 
blocked the whole adaptive process is actually reset. Also it is possible to be more 
generous in the blocking logic - meaning it is possible to allow it to happen more often, 
5 since the filter will fade into a sub optimal, but fairly good, solution and the pitfalls 
normally associated with blocking are avoided. 

4- Alternative approaches to design and implement fixed beam filters 

In parts one and two of this description fixed beamformer 

10 implementations of the invention were presented. In these two implementations one 
simulates a noise structure by placing noise sources in the sphere, then the noise 
covariance matrix is calculated and the optimum filter for that noise structure is 
obtained. In part three of the description a hybrid beamformer implementation was 
discussed. In the third implementation, an adaptive process is employed when there are 

15 significant noises to adapt to, and the fixed solution kicks in when the adaptive process 
in inhibited for some reason (e.g. a strong signal). 

It is proposed here that, assuming one has the infrastructure for an 
adaptive solution, it can be utilized to obtain the fixed solution using the adaptive 
process. For example, lets assume that the adaptive process is implemented on an off 

20 line system using high-level language (like Matlab for example). One can simulate the 
noise structure off-line, i.e. obtain the noise signal on each of the microphones (time 
domain noise sources multiplied by the source steering vector). This noise data can then 
be fed into the simulated off-line adaptive process. Once the adaptive process 
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converges, one can read the final values of the filters' coefficients and use them as the 
optimum solution for the pre-defined noise situation. The disadvantage of using the 
adaptive process in a simulated environment to obtain the fixed weights is that it is time 
consuming. Large data files need to be prepared for the filters to converge and the 

5 adaptive process is a very computation intensive when it is done off-line. Also, the 

existence of an adaptive system simulation has been assumed, and if one does not exist 
it needs to be "prepared.The advantage of this method is that it would provide a more 
accurate solution than the direct methods. The reason is that the direct methods 
determine the gain weight in the frequency domain. It is then necessary to go through a 

10 filter design process that is, by nature, an approximation and includes inherent 

compromises, over which we have no control. Even more so, in the methods discussed 
inparts one and two each filter is designed separately and we have no guaranty that the 
overall beamformer performance (using all the filters concurrently) could not provide a 
better solution. 

15 Running the simulated data through the adaptive process assure us that 

we get the optimum solution for the simulated scenario, that is for the simulated noise 
environment and array structure. For example, if we use the reference channel type 
adaptive filter, the solution will take into account the specific way we actually 
implemented the reference channels - which the separate filter design discussed in part 

20 two does not take into account. 

Another approach proposed in this invention assumes that there is a real 
time working adaptive system. The simulated noise data can be stored on a recording 
media, such as a multi-channel digital tape recorder, or a computer equipped with a 
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multi-channel sound card. The noise data can be injected into the real time working 
system which will converge to the solution, freeze the final filters' coefficients and 
either store them permanently as the fixed solution or transmit them to a hosting system 
to be burned into the fixed beamformer solution. The advantage of this method is that 
5 once the noise data is prepared, the solution is obtained very fast. The adaptive filter 
will converge within seconds. Another advantage of this method is that the fixed 
solution will take into account all kinds of implementation related issues like - fixed 
point and numerical inaccuracies, final dynamic range of the system, differences in the 
input ports of the processor like different AID converters and so on. 

10 Taking the above approach one step further, the present invention 

proposes to create a simulated noise environment using loudspeakers in an an-echoic 
chamber, then running the adaptive system in the chamber and freezing the final values 
of coefficients as the fixed array solution. Loudspeakers are placed in an an-echoic 
chamber to simulate a certain noise scenario - for example two loudspeakers can be 

1 5 placed on each side of the array at 40 degrees and 75 degrees azimuth angle. A 

simulated noise is played through the loudspeakers - for example pink wide band noise. 
The adaptive system runs and converges (within seconds) and the final filters' 
coefficients are stored. The process can be automated - the adaptive system is put in a 
calibration mode, the adaptive system converges and than stores coefficients converged 

20 to as in its own memory as the fixed solution. The calibrated system is than switched off 
from the calibration mode for normal operation. 

The advantage of using the actual working adaptive system is that the 
convergence solution takes into account not just the process itself with all its 
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peculiarities like dynamic range of the processor and the exact implementation of the 
filters, but also unknown factors like the microphones sensitivities and phases, 
mechanical interferences and so on. This is particularly important since it has been 
observed that the fixed solution is very sensitive to some parameters like mismatch in 
5 phases. Also, if the sensors are microphones, for example, and cardioids (uni 

directional) microphones are used instead of omni directional microphones, then the 
mismatch in phase may be such that the actual performance of the filters may be far 
from what was pre-designed. The packaging of the microphone (or other sensor) array 
may also affect the performance strongly. Using the real working adaptive system to 

10 adaptively generate the fixed solution coefficients takes all these parameters into 
account and ensures an optimum solution the given system. 

The disadvantage of the method is that, in general, it is necessary to use 
many simulated noise sources in order to achieve desirable performance improvement. 
Use of one noise source located at one side of an array, for example, may cause the 

15 array to adapt such that the noise source is effectively cancelled while the beam shape 
on the array side opposite is undesirable. However, for a relatively small array, where 
the fixed super directionality is most needed, few noise sources will usually be 
sufficient to provide an improved performance. For instance, in a four cardioids 
microphone array with an aperture of 6" four noise sources are sufficient to provide a 

20 noise rejection of 20 dB at angles over 30 degrees from the look direction. 

An illustrative procedure for generating fixed filter coefficients through 
the use of simulated noise and an actual adaptive system positioned in an an-echoic 
chamber is shown Fig. 3. The first step is to create four random noise files having a 
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white or pink spectrum and a duration of 30 seconds or more (step 74). Next, four 
speakers and an adaptive beamforming system are place in an an-echoic chamber, with 
the angles between the speakers and array look direction being set at -70°, -40°, 40° and 
70° (step 76). The four noise files are fed to the loudspeakers (step 78) and the adaptive 
5 system is allowed to converge to the optimal solution and the filter coefficients 
corresponding to the optimal solution are stored (step 80). 

Another technique to calibrate a system is proposed here. The 
microphone array is placed in the an-echoic chamber and the simulated noise is played 
through the loudspeakers. The output of the array is recorded (no real time DSP system 

10 is present in the chamber). The recorded output is then replayed into the real time 

system. The adaptive process converges and the final filters' coefficients are stored and 
burned into the system as the fixed array solution. This method is sometimes more 
practical when the automatic calibration and burning mechanism is not implemented. It 
is highly inconvenient to perform the down loading and uploading of the coefficient 

15 from a system that is positioned in the chamber. This operation usually requires a 
development system (like In Circuit Emulator or a simulator). It is much more 
convenient to do the recording in the chamber and perform the down loading and 
uploading of coefficients outside were the development system is located. 

An illustrative procedure for generating fixed filter coefficients through 

20 the use of simulated noise, a microphone array positioned in an an-echoic chamber and 
an actual adaptive system positioned outside an an-echoic chamber is shown in Fig. 4. 
As in the procedure of Fig. 3, the first step in the Fig. 4 procedure is to create four 
random noise files having a white or pink spectrum and a duration of 30 (step 82). The 
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next is to place four speakers and an a microphone array in an an-echoic chamber, with 
the angles between the speakers and array look direction being set at -70°, -40°, 40° and. 
70° (step 84). The four noise files are fed to the loudspeakers (step 86) and the 
microphone array's output is recorded on a multi-channel recorder (step 88). The 
5 recorded output is then played into an adaptive beamformer system which is located 
outside the an-echoic chamber and the beamformer is allowed to converge to the 
optimal solution, the coefficients corresponding to the optimal solution being stored for 
use as the fixed filter coefficients (step 90). 

Fig. 5 shows an illustrative procedure for generating simulated noise and 
10 using the simulated noise to generate fixed beamformer coefficients. The first step in 
the procedure is to define the desired noise field scenario and the array configuration 
(step 90). Next, a counting variable indicative of the noise source being considered is 
initialized to one (step 94). A random signal is generated to represent the noise 
emanating from the source under consideration (step 96), and for each sensor, the 
1 5 contribution of noise from the source under consideration is calculated. Calculation of 
noise source contributions involves; initializing to one a counting variable indicative of 
the sensor under consideration (step 98); determining the time delay from the source to 
the sensor under consideration, relative to the time delay to other sensors (step 100); 
and determining the noise source contribution based on the random signal generated in 
20 step 96 and the time delay (step 102). 

After the noise contribution of a source to a particular sensor is 
calculated, a determination is made as to whether all sensors have been considered (step 
104). If all sensors have not been considered, the sensor counting variable is 
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incremented (step 106) and the procedure returns to step 98. When all sensors have 
been considered for a particular source, a determination is made as to whether all . 
sources have been considered (step 108), and if not, the source counting variable is 
incremented (step 1 10) and the procedure returns to step 96. Once the contribution of 
5 each noise source to each sensor has been calculated the generation of the simulated 
noise data is complete. The noise data is then fed to an adaptive procedure which is 
allowed to converge, and the coefficients derived from the converged operation are 
stored for use as the optimal fixed coefficients (step 112). 

1 0 While the present invention has been particularly shown and described in 

conjunction with preferred embodiments thereof, it will be readily appreciated by those 
of ordinary skill in the art that various changes may be made without departing from the 
spirit and scope of the invention. Therefore, it is intended that the appended claims be 
interpreted as including the embodiments described herein as well as all equivalents 

15 thereto. 
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What is claimed is: 

1 . A sensor array for receiving a signal that includes a desired signal and noise, 
comprising: 

a plurality of sensors; 

5 a plurality of filters for filtering the output of each sensor, each filter 

being defined by one or more filter coefficients; and 

a means for combining the outputs of said filters to form a sensor array 

output signal; 

wherein said filter coefficients are adaptively determined so as to 
10 maximize the signal to noise ratio of the array output signal, and wherein upon the 
occurrence of a predetermined event the adaptive determination of said filter 
coefficients is stopped and said filter coefficients are faded into a predetermined set of 
fixed coefficients. 

2. The sensor array as set forth in claim 1, wherein said sensors are 
15 microphones. 

3 The sensor array as set forth in claim 1, wherein said filter coefficients 
are frequency domain coefficients such that for each said filter said frequency domain 
coefficients determine the frequency response of said filter. 

4 The sensor array as set forth in claim 1 , wherein said filter coefficients 
20 are time domain coefficients. 

5 The sensor array as set forth in claim 4, further comprising: 
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a plurality of delay lines, said delay lines corresponding to respective 
outputs of said sensors and receiving respective outputs from said sensors; and 

a plurality of convolvers, corresponding to respective outputs of said 
delay lines, said convolvers being operative to receive respective outputs from said 
5 delay lines and convolve the received outputs with respective filter coefficients to 
generate a plurality of filtered delay line outputs; 

wherein said plurality of filtered delay line outputs are combined by said 
means for combining to form said array output. 

6. The sensor array as set forth in claim 5, further comprising a plurality of 
1 0 signal conditioners for receiving respective outputs from said sensors, sampling the 

received outputs and passing the sampled received outputs to respective delay lines. 

7. The sensor array as set forth in claim 1, wherein said fixed filter 
coefficients are determined by solving directly and non-adaptively an equation 

CT ! v 

w = where C is the noise covariance matrix, v is the steering vector 

°* vCT'v 

15 toward the array look direction, and w op t is a vector having a number of components 
equal to the number of sensors, and where solving said equation for "n" frequencies 
provides "n" coefficients for each filter. 

8. The sensor array as set forth in claim 1 , wherein said fixed filter 
coefficients are determined by placing the array in a simulated noise environment, 

20 letting the adaptation of the filter weights converge to a solution and then setting the 
fixed filter coefficients equal to the coefficients of the converged adaptive solution. 
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9. The sensor array as set forth in claim 1, wherein said fixed filter 
coefficients are determined by simulating a noise environment and the array's response 
to said noise environment, letting the simulated adaptation of the filter weights 
converge to a solution and then storing the coefficients of the converged adaptive 

5 solution for use as the fixed weights of an actual array. 

10. The sensor array as set forth in claim 1, wherein said fixed filter 
coefficients are determined by simulating a noise environment, recording the simulated 
noise generated in said environment, playing back said simulated noise for reception by 
the array, letting the adaptation of the filter weights converge to a solution and then 

10 setting the fixed filter coefficients of the array equal to the coefficients of the converged 
adaptive solution. 

11. An sensor array for receiving signal that includes a desired signal and 
noise, comprising: 

a plurality of sensors; 

15 a delay and sum beamformer for combining the outputs of said sensors to 

generate a beamformer output; 

a reference channel processor for combining the outputs of said sensors 
to generate one or more reference channel signals; 

at least one filter for each said reference channel, each said filter being 
20 defined by one or more coefficients; and 

means for combining the outputs of said filters with said beamformer 
output to generate a sensor array output signal; 
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wherein said reference channel processor and said filters operate to 
maximize the signal to noise ratio of the array output signal, wherein said filter 
coefficients are adaptively determined so as to maximize the signal to noise ratio of the 
array output signal, and wherein upon the occurrence of a predetermined event the 
5 adaptive determination of said filter coefficients is stopped and said filter coefficients 
faded into a predetermined set of fixed coefficients. 

12. The sensor array as set forth in claim 1 1 , wherein said sensors are 
microphones. 

13. The sensor array as set forth in claim 1 1 , wherein said filter coefficients 
10 are frequency domain coefficients such that for each said filter said frequency domain 

coefficients determine the frequency response of said filter. 

14. The sensor array as set forth in claim 1 , wherein said filter coefficients 
are time domain coefficients. 

15. The sensor array as set forth in claim 14, further comprising: 

15 a main channel delay line for delaying the output of said beamformer; 

a plurality of reference channel delay lines, said reference channel delay 
lines corresponding to respective reference channel signals and receiving respective 
reference channel signals; and 

a plurality of convolvers, corresponding to respective outputs of said 
20 reference channel delay lines, said convolvers being operative to receive respective 
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outputs from said reference channel delay lines and convolve the received outputs with 
respective filter coefficients to generate a plurality of filtered delay line outputs; 

wherein said plurality of filtered delay line outputs and said main 
channel delay line output are combined by said means for combining to form said array 
5 output. 

1 6 The sensor array as set forth in claim 1 5, further comprising a plurality 

of signal conditioners for receiving respective outputs from said sensors, sampling the 
received outputs and passing the sampled received outputs to said delay and sum 
beamformer and said reference channel processor. 

10 17. The sensor array as set forth in claim 1 1 , wherein said fixed filter 

coefficients are determined by solving directly and non-adaptively an equation 
Wopt^T'p where C is the noise covariance matrix, p is a vector representing the 
correlation between the output of said beamformer and the output of said reference 
channels, and w op , is a vector having a number of components equal to the number of 

1 5 sensors, and where solving said equation for "n" frequencies provides "n" coefficients 
for each filter. 

1 8. The sensor array as set forth in claim 1 1 , wherein said fixed filter 

coefficients are determined by placing the array in a simulated noise environment, 
letting the adaptation of the filter weights converge to a solution and then setting the 
20 fixed filter coefficients equal to the coefficients of the converged adaptive. 
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19. The sensor array as set forth in claim 11, wherein said fixed filter 
coefficients are determined by simulating a noise environment and the array's response 
to said noise environment, letting the simulated adaptation of the filter weights 
converge to a solution and then storing the coefficients of the converged adaptive 

5 solution for use as the fixed weights of an actual array. 

20. The sensor array as set forth in claim 11, wherein said fixed filter 
coefficients are determined by simulating a noise environment, recording the simulated 
noise generated in said environment, playing back said simulated noise for reception by 
the array, letting the adaptation of the filter weights converge to a solution and then 

10 setting the fixed filter coefficients of the array equal to the coefficients of the converged 
adaptive solution. 

21 . A method of processing a received signal that includes a desired signal 
and noise, comprising the steps of: 

providing an array of sensors; 
15 filtering the output of each sensor through a filter, each filter being 

defined by one or more filter coefficients; and 

combining the outputs of said filters to form a sensor array output signal; 
wherein said filter coefficients are adaptively determined so as to 
maximize the signal to noise ratio of the array output signal, and wherein upon the 
20 occurrence of a predetermined event the adaptive determination of said filter 

coefficients is stopped and said filter coefficients are faded into a predetermined set of 
fixed coefficients. 
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22. The method according to claim 2 1 , wherein said sensors are 
microphones. 

23. The method according to claim 21, wherein said filter coefficients are 
frequency domain coefficients such that for each said filter said frequency domain 

5 coefficients determine the frequency response of said filter. 

24. The method according to claim 21, wherein said filter coefficients are 
time domain coefficients. 

25. The method according to claim 24, further comprising the steps of: 
receiving the outputs of said sensors at a plurality of respective delay 

10 lines; 

receiving the outputs of said delay lines at respective convolvers; 
convolving the received delay line outputs with respective filter 
coefficients to generate a plurality of filtered delay line outputs; and 

combining said plurality of filtered delay line outputs to generate said 

15 array output 

26. The sensor array as set forth in claim 25, further comprising the steps of: 
receiving the outputs of said sensors at respective signal conditioners; 

and 

sampling the received outputs and passing the sampled received outputs 
20 to respective delay lines. 
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27. The sensor array as set forth in claim 21, wherein said fixed filter 
coefficients are determined by solving directly and non-adaptively an equation 

cr'v 

w_ = — where C is the noise co variance matrix, v is the steering vector 

°* vC v 

toward the array look direction, and w opl is a vector having a number of components 
5 equal to the number of sensors, and where solving said equation for * V frequencies 
provides "n" coefficients for each filter. 

28. The sensor array as set forth in claim 21, wherein said fixed filter 
coefficients are determined by placing the array in a simulated noise environment, 
letting the adaptation of the filter weights converge to a solution and then setting the 

10 fixed filter coefficients equal to the coefficients of the converged adaptive solution. 

29. The sensor array as set forth in claim 21, wherein said fixed filter 
coefficients are determined by simulating a noise environment and the array's response 
to said noise environment, letting the simulated adaptation of the filter weights 
converge to a solution and then storing the coefficients of the converged adaptive 

15 solution for use as the fixed weights of an actual array. 

30. The sensor array as set forth in claim 21, wherein said fixed filter 
coefficients are determined by simulating a noise environment, recording the simulated 
noise generated in said environment, playing back said simulated noise for reception by 
the array, letting the adaptation of the filter weights converge to a solution and then 

20 setting the fixed filter coefficients of the array equal to the coefficients of the converged 
adaptive solution. 
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3 J a A method for receiving a signal that includes a desired signal and noise, 

comprising the steps of: 

providing an array of sensors; 

generating a beamformer output by passing the outputs of said sensors 
5 through a delay and sum beamformer; 

generating one or more reference channel signals by passing the outputs 
of said sensors through a reference channel processor; 

filtering each reference channel using at least one filter, each said filter 
being defined by one or more coefficients; and 
10 combining the outputs of said filters with said beamformer output to 

generate a sensor array output signal; 

wherein said reference channel processor and said filters operate to 
maximize the signal to noise ration of the array output signal, wherein said filter 
coefficients are adaptively determined so as to maximize the signal to noise ratio of the 
1 5 array output signal, and wherein upon the occurrence of a predetermined event the 

adaptive determination of said filter coefficients is stopped and said filter coefficients 
are faded into a predetermined set of fixed coefficients. 

32. The method according to claim 3 1 , wherein said sensors are 

microphones. 

20 33. The method according to claim 31, wherein said filter coefficients are 

frequency domain coefficients such that for each said filter said frequency domain 
coefficients determine the frequency response of said filter. 
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34. The method according to claim 3 1 , wherein said filter coefficients are 
time domain coefficients. 

35. The method according to claim 34, further comprising the steps of: 
delaying the output of said beamformer via a main channel delay line; 

5 delaying said reference channel signals via respective reference channel 

delay lines; 

convolving the outputs of said reference channel delay lines with 
respective filter coefficients to generate a plurality of filtered delay line outputs; and 
combining said filtered delay line outputs and said main channel delay 
1 0 line output to generate said array output. 

36. The method according to claim 35, further comprising the steps of: 
receiving the outputs of said sensors at respective signal conditioners; and 

sampling the received outputs and passing the sampled received outputs 
to said delay and sum beamformer and said reference channel processor. 

15 37. The method according to claim 3 1 , wherein said fixed filter coefficients 

are determined by solving directly and non-adaptively an equation Wopt^C^p 
where C is the noise covariance matrix, p is a vector representing the correlation 
between the output of said beamformer and the output of said reference channels, and 
w op t is a vector having a number of components equal to the number of sensors, and 

20 where solving said equation for "n" frequencies provides "n" coefficients for each filter. 
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38, The method according to claim 3 1 , wherein said fixed filter coefficients 

are determined by placing the array in a simulated noise environment, letting the 
adaptation of the filter weights converge to a solution and then setting the fixed filter 
coefficients equal to the coefficients of the converged adaptive. 

5 29. The method according to claim 31, wherein said fixed filter coefficients 

are determined by simulating a noise environment and the array's response to said noise 
environment, letting the simulated adaptation of the filter weights converge to a solution 
and then storing the coefficients of the converged adaptive solution for use as the fixed 
weights of an actual array. 

10 40. The method according to claim 31, wherein said fixed filter coefficients are 

determined by simulating a noise environment, recording the simulated noise generated 
in said environment, playing back said simulated noise for reception by the anay, 
letting the adaptation of the filter weights converge to a solution and then setting the 
fixed filter coefficients of the array equal to the coefficients of the converged adaptive 

15 solution. 

41 . A sensor array for receiving a signal that includes a desired signal and noise, 
comprising: 

a plurality of sensors; 

a plurality of filters for filtering the output of each sensor, each filter 
20 being defined by one or more filter coefficients; and 
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a means for combining the outputs of said filters to form a sensor array 

output signal; 

wherein said filter coefficients are determined by solving an equation 

C~ l v 

w = where C is the noise covariance matrix, v is the steering vector 

°* vcrV 

5 toward the array look direction, and w opt is a vector having a number of components 
equal to the number of sensors, and where solving said equation for "n" frequencies 
provides "n" coefficients for each filter; and 

wherein said noise covariance matrix is determined by defining a spatial 
distribution of noise sources; defining a delay vector for each noise source using said 

1 0 spatial distribution, said delay vector expressing the relative times of arrival of the 

wavefront from said noise source at each sensor; defining a steering vector for each said 
noise source based on said delay vector; using said steering vector to determine the 
contribution of each noise source to said noise covariance matrix; and generating said 
noise covariance matrix by adding the contributions of each noise source and a matrix 

1 5 indicative of spatially distributed white noise. 

42. The sensor array as set forth in claim 41, wherein said equation is solved for **n" 
frequencies to yield "n" coefficients for each filter, and for each filter, said "n" 
coefficients are used to generate filter design values according to an operation selected 
from the group consisting of a Remez Exchange operation and an Inverse Fourier 
20 Transform operation. 
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43. An sensor array for receiving signal that includes a desired signal and 

noise, comprising: 

a plurality of sensors; 

a delay and sum beamformer for combining the outputs of said sensors to 
5 generate a beamformer output; 

a reference channel processor for combining the outputs of said sensors 
to generate one or more reference channel signals; 

at least one filter for each said reference channel, each said filter being 
defined by one or more coefficients; and 
10 means for combining the outputs of said filters with said beamformer 

output to generate a sensor array output signal; 

wherein said filter coefficients are determined by solving an equation 
w op t=C" 1 p where C is the noise covariance matrix, p is a vector representing the 
correlation between the output of said beamformer and the output of said reference 
15 channels, and w opt is a vector having a number of components equal to the number of 
sensors, and where solving said equation for "n" frequencies provides "n" coefficients 
for each filter, and 

wherein said noise covariance matrix is determined by defining a spatial 
distribution of noise sources; defining a delay vector for each noise source using said 
20 spatial distribution, said delay vector expressing the relative times of arrival of the 

wavefront from said noise source at each sensor; defining a steering vector for each said 
noise source based on said delay vector; using said steering vector to determine the 
contribution of each noise source to said noise covariance matrix as measured at the 
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sensors; defining a nulling matrix which indicates how said filter outputs are combined 
to generate said reference channels; determining an array steering vector towards the 
array look direction; determining the contribution of each noise source to each reference 
channel based on said contribution of each noise source at said sensors, said nulling 
5 matrix and said array steering vector; and generating said noise covariance matrix by 
adding the contributions of each noise source to said reference channels and a matrix 
indicative of spatially distributed white noise 

44. The sensor array as set forth in claim 43, wherein said equation is solved 
for 4 V frequencies to yield "n" coefficients for each filter, and for each filter, said "n" 

10 coefficients are used to generate filter design values according to an operation selected 
from the group consisting of a Remez Exchange operation and an Inverse Fourier 
Transform operation. 

45. A method of processing a received signal that includes a desired signal 
and noise, comprising the steps of: 

1 5 providing an array of sensors; 

filtering the output of each sensor through a filter, each filter being 
defined by one or more filter coefficients; and 

combining the outputs of said filters to form a sensor array output signal; 
wherein said filter coefficients are determined by solving an equation 

20 w_ = where C is the noise covariance matrix, v is the steering vector 

vcrV 

toward the array look direction, and Wopt is a vector having a number of components 
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equal to the number of sensors, and where solving said equation for "n" frequencies 
provides "n" coefficients for each filter; and 

wherein said noise covariance matrix is determined by defining a spatial 
distribution of noise sources; defining a delay vector for each noise source using said 

5 spatial distribution, said delay vector expressing the relative times of arrival of the 

wavefront from said noise source at each sensor; defining a steering vector for each said 
noise source based on said delay vector; using said steering vector to determine the 
contribution of each noise source to said noise covariance matrix; and generating said 
noise covariance matrix by adding the contributions of each noise source and a matrix 

1 0 indicative of spatially distributed white noise. 

46 The method according to in claim 45, wherein said equation is solved for 
"n" frequencies to yield "n" coefficients for each filter, and for each filter, said "n" 
coefficients are used to generate filter design values according to an operation selected 
from the group consisting of a Remez Exchange operation and an Inverse Fourier 

15 Transform operation. 

47 A method of processing a received signal that includes a desired signal 
and noise, comprising the steps of: 

providing an array of sensors; 

generating a beamformer output by passing the outputs of said sensors 
20 through a delay and sum beamformer; 

generating one or more reference channel signals by passing the outputs 
of said sensors through a reference channel processor; 
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filtering each reference channel using at least one filter, each said filter 
being defined by one or more coefficients; and 

combining the outputs of said filters with said beamformer output to 
generate a sensor array output signal; 
5 wherein said filter coefficients are determined by solving an equation 

t=c -ip where C is the noise covariance matrix, p is a vector representing the 
correlation between the output of said beamformer and the output of said reference 
channels, and w^t is a vector having a number of components equal to the number of 
sensors, and where solving said equation for "n" frequencies provides "n" coefficients 

10 for each filter; and 

wherein said noise covariance matrix is determined by defining a spatial 
distribution of noise sources; defining a delay vector for each noise source using said 
spatial distribution, said delay vector expressing the relative times of arrival of the 
wavefront from said noise source at each sensor; defining a steering vector for each said 

15 noise source based on said delay vector; using said steering vector to determine the 
contribution of each noise source to said noise covariance matrix as measured at the 
sensors; defining a nulling matrix which indicates how said filter outputs are combined 
to generate said reference channels; determining an array steering vector towards the 
array look direction; determining the contribution of each noise source to each reference 

20 channel based on said contribution of each noise source at said sensors, said nulling 
matrix and said array steering vector; and generating said noise covariance matrix by 
adding the contributions of each noise source to said reference channels and a matrix 
indicative of spatially distributed white noise. 
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48. The method according to claim 47, wherein said equation is solved for 

4 V frequencies to yield "n" coefficients for each filter, and for each filter, said "n" 
coefficients are used to generate filter design values according to an operation selected 
from the group consisting of a Remez Exchange operation and an Inverse Fourier 
5 Transform operation. 
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Define the desired noise fieJd scenario: 

- q directional sources; amplitude s{0 and direction (azimuth and elevation). 

- spatially white noise with amplitude n , 
The array configuration: 

- Sensors location r(i) 

- Null matrix N (The sum of each row of the matrix N should be zero) 

- Look direction (azimuth and elevation). 

The frequency range and resolution to calculate the optimal weights F 
j 



1=1; 
f*F(m) 
0=0; p=0; 
J 



i 

r 



y | =?S i *exp{-ja>T 1 }=3 i exp{ T j27ifRk/c) (the contribution" to the array output) 
1 c is the speed of sound 
k a direction vector pointing the source 

R distance matrix 



C.-XjX,' contribution to the reference channels covarlance matrix. 
x^NdiagMy, 

v is the steering vector towered the look direction 



I 

-JL. 



beam,=v*yi 

p=beam*X 1 correlation vector between the beam out and the references 




Sumation of the contribustion of alt sources 
OlfCj+nl) The reference channel Covariance 
matrix 

p=Z(p | ) The beam to references Co relation 
vector 



So 



si 



1 



w(f)=C- n p The optimal coefficients vector. 




For each reference channel use a filter design program to obtain the time domain 
coefficients that approximate the desired response as defined by W (f) 



7^ 
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Create 4 files of random noise (with a white 
or pink spectrum ) with a duration of at least 
30 second. 



7M 



Place 4 speakers and an adaptive 
beamforming system in anechoic chamber. 
The angles between the speakers direction 
and the array "look direction" is : 
f-70 -40 40 70 ] degrees 



I 



Feed the 4 noise files to the loudspeakers 
(for example by a multichannel sound card, 
or a multichannel tape recorder) 



Let the adaptive system converge to the 
optimal solution and store the optimal 
coefficient in the fixed location. 



fjl — - 
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Create 4 files of random noise (with a white 
or pink spectrum ) with a duration of at least 
30 second. 



i 



P lace 4 soeakers and a *n«creph*n^ 

a^rc^Kf in QQ^choic chamber. 
The angles between the speakers direction 
and the array "took direction" is : 
f-70 -40 40 70 1 degrees 



Feed the 4 noise files to the loudspeakers 
(for example by a multichannel sound card, 
or a multichannel tape recorder) 



62- 



R^rJ) microphone array 

output on a multi- 
channel recorder. 



Feed the recorded data to a real time system 
that implements the adaptive beamformer * 
and let it converge. 

Capture the coefficients and store them as 
the optimal fixed coefficients. 



4 
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Define the desired noise field scenario: 

. q directional sources with amplitude s(i) and direction fa® f(i)] (azimuth and elevation). 

- spatially white noise with amplitude n, 
The array configuration: 

- Sensors location r(i) 

- NuU matrix N (The sum of each row of the matrix N should be zero) 

- LooK direction (azimuth and elevation). 



r=1; (source l) 



MY 



Generate a random signal a t (t) 



(sensor J) 



92 



t^-r, )k/c (time delay) 
k a direction vector pointing the source 

|sin(e)cos«>),sin(e)sin(4») f cos(e)r 
r sensors location vector 



y 0 (t)^n(t)+s,(t-V 

Accumulate the cor 



contribution of that source I to each element j 



M*1 



r=M k 




9£ 



700 



Feed the simulated data to an adaptive procedure that implementes the adaptive 
beamformer and let it converge 

Store the final results as the optimal fix coefficients 



hill. 
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